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Abstract 

We  propose  a  Hermite  representation  for  deformable  contour  finding.  This  representa¬ 
tion  compares  favorably  in  terms  of  versatility  and  controllability  with  other  local  contour 
representations  that  have  been  used  previously  for  this  purpose.  The  Hermite  representation 
allows  a  compact  representation  of  curved  shapes,  without  the  smoothing  out  of  corners.  It 
is  also  well  suited  for  both  interactive  and  tracking  applications. 

The  Hermite  representation  is  used  to  formulate  the  contour  finding  problem  as  an  opti¬ 
mization  problem  using  a  maximum  a  posteriori  energy  criterion.  Optimization  is  performed 
by  dynamic  programming.  Our  approach  to  contour  tracking  decouples  the  effects  of  trans¬ 
formation  and  deformation,  using  a  template  matching  strategy  to  robustly  account  for  the 
transformation  effect.  We  demonstrate  these  ideas  on  a  variety  of  images  from  different 
domains. 
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1  Introduction 


Image  segmentation  by  boundary  finding  is  one  of  the  central  problems  in  computer  vi¬ 
sion.  This  is  because,  among  the  features  that  can  be  used  to  distinguish  objects  from  their 
backgrounds,  such  as  color  and  texture,  shape  is  usually  the  most  powerful.  For  detect¬ 
ing  instances  of  objects  with  fixed  and  known  shape,  the  generalized  Hough  transform  or 
template  matching  is  well  suited  (see  [3]).  For  cases  where  there  exists  some  flexibility  in 
the  object’s  shape  (either  w.r.t.  a  previous  frame  in  a  tracking  application,  or  w.r.t.  a  user- 
supplied  shape  in  an  interactive  object  delineation  setting)  deformable  contour  models  have 
found  widespread  use. 

Deformable  contours  (also  called  active  contour  models,  or  “snakes”)  are  energy¬ 
minimizing  models  for  which  the  minima  represent  solutions  to  contour  segmentation  prob¬ 
lems.  They  can  overcome  problems  with  traditional  bottom-up  segmentation  methods,  such 
as  edge  gaps  and  spurious  edges,  by  the  use  of  an  energy  function  that  contains  shape  infor¬ 
mation  in  addition  to  terms  determined  by  image  features.  The  additional  shape  information 
can  be  seen  as  a  regularization  term  in  the  fitting  process.  Once  placed  in  image  space,  the 
snake  deforms  to  find  the  most  salient  contour  in  its  neighborhood,  under  the  influence  of 
the  generated  potential  field. 

Extensive  work  heis  been  reported  on  deformable  contours  since  their  emergence  in  the 
late  eighties;  examples  are  [1,  4-6,  8-12,  14-16].  A  useful  way  to  characterize  the  different 
approaches  is  along  the  following  dimensions: 

•  contour  representation 

•  energy  formulation  (internal  and  external) 

•  contour  propagation  mechanism  (spatial  and  temporal) 

In  Section  2  we  review  the  various  contour  representations  that  have  been  used.  We 
proposed  a  new  local  representation  for  the  deformable  contour  framework,  based  on  Hermite 
interpolating  cubics;  see  Section  3.  As  will  become  apparent,  it  has  several  advantages;  in 
particular,  it  handles  both  smooth  and  polygonal  curves  naturally. 
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We  formulate  the  solution  to  the  contour  finding  problem  by  a  maximum  a  posteri¬ 
ori  (MAP)  criterion.  This  leads  to  an  internal  energy  formulation  which  contains  squared 
terms  of  deviations  from  the  expected  Hermite  parameter  values.  The  external  energy  terms 
describe  the  typical  image  gradient  correlations;  see  Section  4.1.  The  resulting  energy  mini¬ 
mization  is  performed  by  dynamic  programming,  which  gives  the  optimal  solution  to  contour 
finding  for  a  certain  search  region;  see  Section  4.2. 

One  of  the  well-known  limitations  of  deformable  contours  is  that  their  initial  placement 
has  to  be  close  to  the  desired  object  boundary  in  order  to  converge.  In  tracking  applications, 
this  assumption  might  be  violated.  To  keep  the  problem  computationally  tractable,  we 
propose  to  decouple  the  effects  of  transformation  and  deformation;  see  Section  4.3. 

Experiments  on  a  variety  of  images  are  presented  in  Section  5,  and  conclusions  in  Sec¬ 
tion  6. 

2  Contour  Representations — A  Review 

Contour  representations  can  be  roughly  divided  into  two  classes,  depending  on  whether  they 
are  global  or  local.  Global  representations  are  those  where  changes  in  one  shape  parameter 
affect  the  entire  contour,  and  conversely,  local  change  of  the  contour  shape  affects  all  pa¬ 
rameters.  These  representations  are  typically  compact,  describing  shape  in  terms  of  only  a 
few  parameters.  This  is  an  advantage  in  a  recognition  context,  i.e.  when  trying  to  recover 
these  parameters  from  images,  because  of  lower  complexity.  A  useful  class  of  shapes  easily 
modeled  by  a  few  global  parameters  axe  the  super-quadrics  [15],  which  are  generalizations 
of  ellipses  that  include  a  degree  of  “squareness”.  To  these  shapes,  one  can  add  global  defor¬ 
mations,  such  as  tapering,  twisting  and  bending  [2].  A  more  general  global  representation  is 
the  Fourier  representation  [14].  It  expresses  a  parameterized  contour  in  terms  of  a  number 
of  orthonormal  (sinusoidal)  basis  functions.  Arbitrary  contours  can  be  represented  in  any 
detail  desired,  given  a  sufi&cient  number  of  basis  functions. 

Local  representations  control  shape  locally  by  various  parameters.  This  flexibility  makes 
local  representations  well  suited  in  a  shape  reconstruction  context,  as  is  the  case  when  de¬ 
forming  a  contour  to  fit  image  data.  The  simplest  contour  representation  is  an  ordered  list  of 
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data  points.  More  compact  representations  describe  contours  in  terms  of  piecewise  polyno¬ 
mials.  Each  segment  of  the  parameterized  contour  (xj(i),  yi{t))  is  described  by  a  polynomial 
in  t.  The  lowest-degree  interpolating  polynomial  is  of  degree  one,  leading  to  a  contour  rep¬ 
resentation  by  polylines  and  polygons.  More  flexibility  is  possible  by  the  use  of  higher  order 
polynomials,  generally  cubic  polynomials;  they  are  the  lowest  degree  polynomials  for  which 
derivatives  at  the  endpoints  can  be  specified.  Higher-order  polynomials  tend  to  bounce  back 
and  forth  in  a  less  controllable  fcishion  and  therefore  are  used  less  frequently  for  interpolation 
purposes. 

Natural  cubic  splines  are  piecewise  third-degree  polynomials  which  interpolate  control 
points  with  (7°,  and  C'^  continuity.  The  natural  cubic  spline  parameters  depend  on 
all  the  control  points,  which  makes  them  a  global  representation.  B-splines,  on  the  other 
hand,  have  a  local  representation,  where  contour  segments  depend  only  on  a  few  neighboring 
control  points.  This  comes  at  the  price  of  not  interpolating  the  control  points.  The  same  C®, 

and  continuity  as  with  natural  splines  is  now  achieved  at  the  join  points  of  connecting 
segments.  By  replicating  control  points,  one  can  force  a  B-spline  to  interpolate  the  control 
points.  A  last  interesting  property  is  that  the  B-spline  can  be  specified  so  that  it  performs 
a  least-squares  fit  on  the  available  data  points. 

In  previous  work,  three  local  representations  have  been  used  for  deformable  contour  find¬ 
ing:  point  representations,  polygonal  chains,  and  uniform  B-splines.  These  representations 
have  the  following  disadvantages  when  used  for  the  contour  finding  task. 

Manipulating  contours  on  the  fine  scale  offered  by  pixel-by-pixel  representations  leads 
typically  to  high  computational  cost  (for  example,  note  the  high  complexity  incurred  in  [8]). 
The  incorporation  of  a-priori  shape  information  in  this  featureless  representation  is  difficult. 
If,  on  the  other  hand,  a  contour  is  represented  by  a  few  (feature)  points,  and  contour  finding 
considers  only  image  data  in  the  local  neighborhoods  of  these  points,  no  use  is  made  of  data 
at  intermediate  locations,  which  makes  the  approach  prone  to  image  noise. 

The  polygonal  chain  representation  [6]  overcomes  some  of  these  problems.  However,  it  is 
not  well  suited  to  representing  curved  objects,  requiring  many  control  points  to  be  adequate. 
In  an  interactive  object  delineation  setting,  this  is  tedious.  For  tracking  applications,  the 
placement  of  control  points  close  to  each  other,  typical  also  of  point  representations,  leads  to 
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stability  problems.  This  is  because  for  most  contour  finding  approaches  using  local  represen¬ 
tations,  a-priori  shape  information  is  encoded  for  each  control  point  w.r.t.  to  its  neighboring 
control  points  (i.e.  curvature  [9,  12,  16],  affine  coordinates  [10]).  If  control  points  are  close 
together,  small  deviations  due  to  image  noise  or  contour  propagation  will  result  in  large 
changes  of  local  shape  properties. 

B-splines  present  an  efficient  and  natural  way  to  represent  smoothly  curved  objects.  For 
objects  with  sharp  corners  they  are  less  well  suited;  the  continuity  smooths  out  any  regions 
of  high  curvature  on  a  contour.  The  fact  that  B-splines  do  not  interpolate  the  control  points 
can  be  considered  a  drawback  in  an  interactive  object  delineation  setting  (think  of  a  physician 
pointing  to  specific  locations  in  medical  images).  The  above-mentioned  use  of  control  point 
duplication  can  take  care  of  this,  but  then  straight  line  segments  appear  around  the  new 
(7°  continuous  control  points.  Without  user  intervention,  when  to  duplicate  control  points 
becomes  a  difficult  decision;  for  example,  Menet  Jll]  duplicates  control  points  in  regions 
where  after  M  steps  of  contour  deformation,  the  curvature  is  higher  than  a  user-supplied 
threshold  0. 


3  The  Hermite  Representation 


The  previous  considerations  lead  us  to  propose  the  Hermite  representation  for  deformable 
contour  finding.  Hermite  contours  are  piecewise  cubic  polynomials,  which  interpolate  the 
control  points  po,...,PN-  In  each  interval,  the  Hermite  cubic  Q{s,t)  =  [a:(s,t)  ^(s,t)]  is 
specified  by  the  positions  pi_i,  pi  and  tangent  vectors  t-~  at  the  endpoints. 

Let  Q  be  an  arbitrary  cubic  polynomial 

Q  =  T-C  (1) 


where 


T  =  [t^  1] 


C  = 


0/x 

hx  by 

Cx  Cy 

dx  dy 
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with  tangent  vector  Q'(t) 


Q'  =  T'  •  C  =  [3t^  2t  1  0]  •  C 


(2) 


Given  the  Hermite  parameter  matrix 

Hi  =  [hi^hiy]  =  [pi_i  Pi  Tjti  T.-f  (3) 

the  corresponding  Hermite  coefficient  matrix  Ch-  can  be  derived  [7]  as 


We  collect  all  the  Hermite  parameters  in  state  vector  H  for  later  use: 


When  considering  the  same  criteria  of  usefulness  for  the  contour  finding  problem  as  we 
discussed  in  the  previous  section  for  the  point-,  polygon-  and  spline-based  representations, 
we  note  that  the  Hermite  representation 

•  can  efficiently  represent  both  smooth  and  sharp  contours.  This  is  because  smooth  con¬ 
tours  are  well  represented  by  the  Hermite  interpolating  cubics,  while  at  the  same  time, 
arbitrarily  sharp  corners  can  be  easily  generated  at  the  control  points  by  adjustment 
of  the  left  and  right  tangent  vector  parameters 

•  interpolates  the  control  points 

•  is  explicit  in  those  features  that  can  be  measured  from  image  data:  position  and 
direction  of  gradient  at  control  points.  This  allows  us  to  prune  the  search  space  during 
contour  finding,  as  we  will  see  in  the  next  section. 
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4  Contour  Detection 


4.1  MAP  formulation 


A  maximum  a  posteriori  (MAP)  criterion  is  formulated  for  the  solution  of  the  contour  finding 
problem.  The  aim  is  to  find,  from  among  all  possible  contours,  the  contour  that  matches  the 
image  data  best  in  a  probabilistic  sense.  Let  d(a:,j/)  be  the  original  normalized  image  and 
y)  the  image  template  corresponding  to  the  Hermite  parameters  H.  We  want  to  find 
Hmap  which  maximizes  the  probability  that  tn  occurs  given  d,  e.g.  „  |d).  tui., .  t> 

is  then  the  maximum  a  posteriori  solution  to  the  problem.  Bayes’  rule  gives 


=  m^xP(tHld) 


H 


=  max 

H 


P(d|tH)  P(tH) 
P(d) 


(5) 


where  P(d|tH)  is  the  conditional  probability  of  the  image  given  the  template,  and  P(tH) 
and  P(d)  are  the  prior  probabilities  for  the  template  and  image,  respectively.  Taking  the 
natural  logarithm  on  both  sides  of  eq.  (5)  and  discounting  P(d),  which  does  not  depend  on 
H,  leads  to  an  equivalent  problem  of  maximizing  objective  function  U : 


=  nmx(lnP(tH) +  lnP(d|tH)) 


H 


This  equation  describes  the  trade-off  between  a-priori  and  image-derived  information. 

If  the  image  is  regarded  as  a  noise-corrupted  template  with  additive,  independent  noise 
that  is  zero-mean  Gaussian,  we  have  P(d|tH)  =  P(d|tH  -|-  n)  =  P(n|d  —  tu);  thus 


and 


^(ditH)  =  n 

V  27rcr„ 


In  P(d|tH)  =  constant  +  ^ 

tjjla:,!/) 


(‘i(x,y)-tn(x,y)Y 


{d{x,y)-tH{x,y)y 

2a? 


(6) 


(7) 


The  last  term  can  be  replaced  by  a  correlation  term  d-tn,  approximating  ||d||  and  lltnll 
by  constants.  For  ||d||  =  1  and  lltuH  =  1  we  obtain 


max  In  P(d|tH) 
H 


min  (1  -  d  •  tn)  =  rnin  Pext 

XI  H 


(8) 
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A  similar  derivation  was  described  by  Rosenfeld  and  Kak  [13]  and  Staib  and  Duncan  [14]. 
We  model  the  prior  probability  for  a  Hermite  contour  H  as 

P(H)  =  P(H|H)  =  constant  •  n '  ^  (9) 

i 

where  H  represents  an  expected  contour.  H  is  typically  obtained  as  the  sample  mean  of 
contours  generated  in  a  training  phase,  or  as  the  contour  obtained  by  prediction  during 
tracking.  Ci  acts  as  a  weighting  measure  for  the  various  dimensions.  In  case  of  an  open 
contour,  we  set  the  cr’s  of  Tq  and  to  infinity. 

In  tracking  applications  the  contour  typically  undergoes  a  tranformation  T  (for  example, 
translation,  rotation  and  scale)  for  which  one  does  not  want  to  penalize.  Our  modeling  pro¬ 
cess  assumes  that  any  transformation  of  the  contour  which  one  does  not  want  to  penalize  has 
already  been  performed  before  eq.  (9)  is  applied.  Any  further  contour  change  is  considered 
as  a  deformation  from  an  expected  contour  and  is  thus  penalized. 

Taking  the  natural  logarithm  gives 

max  InP(tH)  =  min  =  min  Pint  (10) 

H  ^  ^  H  “  2af  H  ^  ^ 

4.2  Dynamic  Programming 

There  are  many  ways  to  solve  the  minimization  problem 

imn  E  =  imn  (Pint  -f  Pext).  (11) 

H  H 

Variational  calculus  methods  have  been  used  extensively  for  continuous  parameter  spaces 
where  derivative  information  is  available  [5,  9,  11,  12,  14,  15].  For  discrete  search  spaces 
one  possibility  is  to  use  A.I.  search  techniques  (e.g.  best-first,  simulated  annealing,  genetic 
algorithms).  We  will  use  a  discrete  enumeration  technique  based  on  dynamic  programming 
(DP)  which  was  popularized  by  Amini  et  al.  [1],  and  used  since  in  [8, 10].  The  advantages  of 
dynamic  programming  w.r.t.  variational  calculus  methods  are  in  terms  of  stability,  optimality 
and  the  possibility  of  enforcing  hard  constraints  [1].  For  dynamic  programming  to  be  efficient 
compared  to  the  exhaustive  enumeration  of  the  possible  solutions,  the  decision  process  should 
be  Markovian.  This  is  typically  the  case  if  the  a-priori  shape  component  Pint  contains  a 
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summation  of  terms  which  depend  only  on  parameters  that  can  be  derived  locally  along  the 
contour. 

For  the  case  of  open  contours,  our  objective  function  can  be  written  as 

E  =  F;i(po,r^,ri  ,pi)  +  •••  + 

•S'Ar(PN-li Pn)  (12) 

Applying  the  dynamic  programming  technique  to  our  formulation  involves  generating 
a  sequence  of  functions  of  two  variables,  Si  with  i  =  0, ...  ,A^  —  1,  where  for  each  Si  a 
minimization  is  performed  in  two  dimensions.  Si  are  the  optimal  value  functions  defined  by 

5o(ti”,Pi)  =  min  £;i(po,T^,ri",pi) 

Po-^o 

5i(7-i;i,Pi+i)  =  min(si_i(pi,r/^)  +  Ei(pi,r^,Tr;.i,pi+i)) 

Pi>Y 

i  =  1,...,  A^-1  (13) 

If  Pi  and  Tj~  (rA )  range  over  Np  and  Np  values  at  each  index  i,  the  complexity  of  the 
proposed  algorithm  is  0{NNpNj). 

The  above  formulation  is  for  open  contours.  For  closed  contours,  where  the  first  and  last 
control  points  are  equal,  we  apply  the  same  algorithm  as  for  the  open  contour  case,  and 
repeat  it  for  all  Np  possible  locations  of  the  first  (last)  control  point,  while  keeping  track  of 
the  best  solution.  The  complexity  increases  to  0{NNpN^). 

Speedup  can  be  achieved  by  a  multi-scale  approach.  Here  contour  finding  is  first  done 
on  a  lower  resolution  image  to  find  an  approximated  contour.  This  can  be  done  with  a 
coarse  discretization  of  the  parameter  space  (i.e.  requiring  smaller  Np  and  Np  for  the  same 
parameter  range).  At  the  finer  level,  the  originally  desired  discretization  can  be  achieved  by 
decreasing  the  parameter  range  to  lie  around  the  solution  found  at  the  coarse  level. 

At  each  scale,  the  algorithm  can  be  speeded  up  by  discounting  improbable  control  point 
locations  before  starting  the  DP  search.  A  measure  of  “improbability”  can  be  specified 
in  terms  of  weak  image  gradient  strength  or  the  dot  product  between  the  measured  and 
expected  gradient  directions  (the  latter  are  explicit  in  the  Hermite  representation).  If  all 
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the  candidate  control  point  locations  are  rated  similarly  (e.g.  standard  deviation  of  ratings 
below  a  threshold),  it  is  more  robust  to  consider  all  of  them. 

T  For  closed  contours,  one  can  use  only  single-pass  DP  for  the  closed  contour  and  optimize 

for  the  remaining  Eo(pn>7-N’''’o  >Po)  while  assigning  to  po  and  PN  the  optimal  values  found 
for  the  open  contour  case.  Of  course,  all  these  speedup  procedures  lose  the  optimality 
property  of  DP.  Nevertheless,  the  last  two  methods,  which  we  implemented,  performed 
satisfactorily  in  practice. 

4.3  Contour  Tracking 

The  high  computational  cost  of  dynamic  programming,  and  of  other  search  methods  which 
do  not  get  stuck  at  the  closest  loccil  minimum,  makes  search  feasible  only  in  a  limited 
neighborhood.  For  interactive  contour  delineation  this  is  fine,  since  the  user  is  likely  to  place 
well-positioned  control  points,  very  close  to  the  desired  contour.  In  tracking  applications 
this  requirement  is  often  unrealistic.  On  the  other  hand,  it  is  our  observation  that  the  effects 
of  deformation  are  often  small  from  frame  to  frame  once  rigid  motion  is  accounted  for. 

We  therefore  decouple  the  effects  of  motion  and  deformation  on  the  contour,  searching 
first  for  transformation  parameters  T  =  [t,  s]  with  t,  <!>  and  s  denoting  translation,  rotation 
and  scaling.  T  is  found  w.r.t.  the  undeformed  contour,  after  which  search  continues  for  the 
deformation  parameters.  The  first  stage  is  robustly  performed  by  template  matching  (or  a 
Hough  transform  [10])  on  a  Gaussian-blurred  gradient  image.  The  second  stage  is  the  DP 
approach  described  earlier.  Both  stages  use  motion  prediction  methods;  template  matching 
at  time  t-M  searches  in  a  parameter  range  centered  around  predicted  transformation  r(t-|-l) 
using  predicted  template  H(t  -f  1).  H(t  -t- 1)  is  also  the  initial  contour  of  the  DP  search. 

For  simplicity,  we  currently  use  T{t  -f  1)  =  T{t)  and  H(t  +  1)  =  H(t).  More  generally, 
T{t  -f  1)  =  p{t  +  1)  where  p  is  a  best-fitting  nth  order  polynomial  at  {t  —  M,  T{t  —  M)), . . . , 

4  {t  —  l,t).  Similar  considerations  hold  for  H(t  +  1).  If  the  time  span  M  over  which  an  nth 

order  model  holds  is  large,  it  is  efficient  to  use  a  recursive  predictor  such  as  a  Kalman  filter. 
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5  Experiments 


We  have  performed  experiments  with  the  proposed  combination  of  Hermite  representation 
and  template-plus-DP  search  in  both  interactive  and  tracking  settings.  The  associated  tem¬ 
plate  matching  parameters  were  range  and  discretization  of  the  transformation  parameters 
(translation,  rotation  and  scale).  DP-related  parameters  included  the  initial  values  of  the 
Hermite  parameters,  their  range  and  discretization,  as  well  as  the  weighing  parameters.  The 
locations  considered  around  control  point  pi  lay  on  a  rectangular  grid  with  x-axis  perpen¬ 
dicular  to  pi^.1  —  Pi_i.  The  Hermite  gradients  rj  were  described  in  terms  of  length  and 
direction  4>i.  Typically,  N  =  5,  Np  =  9  (iVp  =  4  after  pruning),  iV)  =  3,  =  9. 

Figure  1  demonstrates  the  versatility  of  the  Hermite  representation.  Different  initial 
contours  are  placed  by  the  user  as  shown  in  Figure  la.  Figure  lb  shows  the  search  region 
covered  by  DP  for  the  initial  control  point  placement;  for  each  contour  segment  the  Hermite 
cubics  are  shown  corresponding  to  and  <f>i+irr,iJ  for  fixed  (initial) 

control  point  locations  and  I  =  /max-  Many  different  Hermite  contours  which  lie  within 
this  search  region  are  not  displayed.  Figure  Ic  shows  the  result  of  contour  finding  by  DP. 
One  observes  a  wide  variety  of  shapes  that  have  been  accurately  described  by  the  Hermite 
representations,  from  the  smoothly  varying  contour  of  the  mug  rim  to  the  sharp  corners  of 
the  square  pattern,  with  a  curved  horizontal  segment  joining  at  the  corner.  They  compare 
favorably  with  possible  representations  by  polygonal  chains,  splines,  or  Fourier  descriptors. 

For  completeness,  we  also  show  in  Figure  Id  the  conditioned  Sobel  gradient  image,  which 
is  used  by  the  DP  algorithm.  We  use  a  conditioned  image  instead  of  the  original  Sobel  image 
in  order  to  amplify  weak  but  probable  edges.  This  is  done  based  on  local  considerations, 
taking  into  account  the  mean  fj,  and  standard  deviation  cr  in  an  n  x  n  neighborhood.  A 
linear  remapping  is  applied  to  the  image  data  at  {x,y)  if  a  is  greater  than  a  user-specified 
threshold. 

Figure  2  shows  different  instances  of  initial  placement  and  contour  detection  on  an  MR 
image  of  the  human  brain.  Finally,  Figure  3  shows  a  tracking  sequence  of  a  head  using  the 
proposed  combination  of  coarse-scale  template  matching  and  DP. 
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(C)  (d) 


Figure  1;  Mug  image:  (a)  contour  initialization,  (b)  search  region,  (c)  contour  detection,  (d) 
conditioned  Sobel  image 

6  Conclusion 


We  have  proposed  the  use  of  a  Hermite  representation  for  deformable  contour  finding.  We 
have  shown  it  to  have  advantages  over  point-,  polygonal-  and  spline-based  representations  in 
terms  of  versatility,  stability  and  controllability.  A  decoupled  approach  to  contour  tracking 
was  proposed  based  on  template  matching  on  a  coarse  scale  to  account  for  motion  effects, 
and  a  DP  formulation  on  a  finer  scale  to  account  for  deformation  effects. 
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Figure  2:  Brain  MR  image:  (a)  contour  initialization,  (b)  contour  detection 


Figure  3:  A  tracking  sequence  of  a  head  {t  =  0, 8, 24, 30) 
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